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Abstract We review our quasiparticle model for the thermodynamics of strongly interacting 
matter at high temperature, and its extrapolation to non-zero chemical potential. 
Some implications of the resulting soft equation of state of quark matter at low 
temperatures are pointed out. 

Keywords: QCD equation of state, quasiparticles, quark stars 



1. Introduction 

The question for the equation of state of strongly interacting matter is a 
link between many-particle physics and astrophysics/cosmology. Calculated 
by means of statistical quantum field theory, it serves as a necessary input, e. g., 
in models of the early universe, or in the context considered here, it determines 
the structure of stars. Astrophysical observations, such as the mass and the 
radius of dense stars, may in turn also impose constraints on the equation of 
state of deconfined quark matter. 
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We will focus on static spherically symmetric stars which are described by 
the Tolman-Oppenheimer-VoUcov equations. At low densities, up to a few times 
nuclear density no, matter consists of interacting hadrons. Theoretical models 
for this state have to start from various assumptions, as for the included states 
and their interactions. Naturally, the results for the hadronic equation of state 
become notably model dependent at densities exceeding approximately 2no. 
This is reflected in uncertainties of the predictions for the shell structure of neu- 
tron stars, cf. [1]. At some higher energy density, hadronic matter undergoes a 
phase transition or a crossover to a quark-gluon plasma (at high temperatures), 
or a color-superconducting state (at low temperatures). Then the system can 
be described directly in terms of the fundamental degrees of freedom - quarks 
and gluons. Notwithstanding, the coupling strength is still large in the regime 
of physical interest, and perturbative QCD is not rehable or, at least, the cal- 
culations have to be interpreted very carefully. Calculations based on various 
effective theories, on the other hand, rest again on assumptions which seem 
hard to control a priori. 

Non-perturbative results have been obtained from first principles by lattice 
QCD computations at zero chemical potential and temperatures up to a few 
times the transition temperature Tc. At low temperatures and for non-zero 
chemical potential ^, as relevant for dense stars, Monte-Carlo calculations are, 
however, hampered by the sign problem. Several approaches to cope with it 
have been proposed only recently. The available results still have rather large 
uncertainties, and they do not yet cover the range of temperatures and chemical 
potential required in the present context. 

This makes worthwhile an approach which extrapolates, with as few assump- 
tions as possible, lattice QCD data from zero chemical potential to /x > 0. In 
the following we will outline a thermodynamical quasiparticle model, which 
can be derived in a series of approximations. 

2. Resummation and quasiparticle models 

2.1 0^ theory 

For the sake of simplicity we consider, for the moment being, a macroscopic 
system described by a scalar field theory. Following [2], the thermodynamic 
potential, Q = ~pV, can be calculated from a functional of the full propagator 

A = (Ao-i-n)-\ 

^ = \f (in(-A-i) + ha) - $[A] , (1) 

evaluated at the stationary point 50./ 5 IS. = 0. This corresponds to a self- 
consistent calculation of the self-energy, 11 = 2(5$ / (5A, where $ is the series 
of 2-particle irreducible diagrams. 
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In this scheme, thermodynamically self-consistent approximations [3] can 
be derived by truncating the expansion of which amounts to a resummation 
of whole classes of diagrams in perturbation theory. To leading order, for the 
massless case and in the MS scheme in 4 — 2€ dimensions. 
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where nb{x) = (e^ — and a; = (A;^ + H)^^^, the self-energy is simply a 
mass term. The non-perturbative gap equation requires the renormaUzation of 
the bare coupling qq. The resummation of the set of 'chain' graphs yields for 
the coupling at the scale of Mandelstam s 
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Expressing now 50 by g{s) leads to the well-defined relation 
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whose solution is interpreted as a temperature dependent quasiparticle mass 
squared. The corresponding pressure reads 
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For both calculational details as well as for a discussion of this approximation 

and its relation to other approaches we refer to [4]; here we only emphasize 
its structure. The first term on the rhs is simply the pressure of free massive 
particles. Written in the form p = p"^{T, m) — B(T), the function B is related 
to m{T) such that the entropy s = dp/dT reduces, due to the stationarity of 
ri, to the entropy s'^(r, m) of an ideal gas. 

2.2 QCD at finite temperature 

HTL quasiparticle model. In QCD, the truncation of a resummation scheme 
based on 2-point functions is delicate because of gauge invariance. However, it 
can be argued that appropriate approximate solutions of the Schwinger-Dyson 
equations can yield reasonable approximations for O, or the pressure, see the 
schematic Figure 1 . Putting this to a test, we consider the representative case 
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■^exact -^approx ^ 

Figure 1. The exact and the truncated functional for the thermodynamic potential. For the 
latter, a propagator 'near' to the stationary point (which itself may be unphysical) can give a 
physically reasonable approximation. 



of the pure gauge plasma, 
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where the traces are taken over the group indices and the 4-momentum, and 
S = Gq^ — G^^ is the ghost self-energy. The propagators can be approximated 
by the hard thermal loop (HTL) contributions, which are gauge invariant and 
have the correct limit for hard momenta near the light cone. The resulting 
approximation for the pressure [5], in condensed form 



Va 



Tr 



In(-Z), 



+ TY[ln(-Go^)], 



(7) 



has similar properties as the corresponding expression in the scalar theory. 
With the standard 2-loop running coupling, it agrees with the lattice data for 
temperatures down to 3Tc, see Figure 2, which is a noteworthy improvement 
compared to the perturbative results. Similar results have been obtained by 
calculating directly the HTL-resummed entropy [7]. 



Phenomenological quasiparticle model. Taking into account only the dom- 
inant contributions in (7), namely the quasiparticle contributions of the trans- 
verse gluons as well as the quark particle-excitations for 7^ 0, we arrive at 
the quasiparticle model [8]. The dispersion relations can be even further sim- 
plified by their form at hard momenta, ujf = k"^ + ml, where nii ~ gT are the 
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Figure 2. Comparison of the pressure, in units of tfie free pressure, of the SU(3) plasma from 
the HTL quasiparticle approximation (7) vs. lattice data [6]. 



asymptotic masses. With this approximation of the self-energies, the pressure 
reads in analogy to the scalar case 



piT) = J2p"'{T,m,)-BiT). 



(8) 



Conceding an enhancement of the running coupling in the infrared, parameter- 
ized by Ts > in an ansatz compatible with the pertuibative limit. 



2 ' 



(9) 



the thermodynamic lattice data can be quantitatively described even down to T^, 
for an example see Figure 3. The coupling g{T) obtained from available lattice 
QCD data with a given number of flavors will be an input for the extrapolation 
to non-zero chemical potential, as outlined in the following. 

2.3 Non-zero chemical potential 

The phenomenological quasiparticle model can be generalized to non-zero 
chemical potential, where the quasiparticle masses of the gluons and quarks 
read 
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TIT, 

Figure 3. The quasiparticle fit (solid line) of the lattice data [9] (open symbols) for the pressure 
in QCD with Nf — 2 light flavors. The full symbols, representing data with large quark masses, 
agree with the results for the pure SU(3) plasma (hatched band); for details see [10]. 
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and the pressure, analogous to Equation (8), now depends also on /i. The 
predictive power of the approach becomes obvious by noting that the depen- 
dance of the coupling on T and jjL is completely governed by the requirement 
of thermodynamic consistency [11]: Maxwell's relation, 

ds{T,fj.,mi) _ dn{T,fi,mi) 
~ dT ' 

where n = dp/d^is the particle density, implies 



dT 



(11) 



While nonlinear in (the coefficients a^^x and b are lengthy integral expres- 
sions), the partial differential equation is linear in the derivatives. It can thus 
be solved by the method of characteristics, with the boundary conditions given 
by the coupling at /i = 0, as obtained from finite-T lattice QCD. The resulting 
elliptic flow shown in Figure 4, which maps the equation of state from /i = 
to > 0, is plausible from the physical intuition. In the perturbative limit, 
T cfj,/-K with c = ( ^^9jv^^^ ^ 1. It is noteworthy that the correspon- 
dence r ~ /i/TT holds with a good accuracy even when g is not small. A similar 
observation was made in the HTL quasiparticle approach [12]. 
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Figure4. The characteristics of the flow equation (1 1) for A'^/- = 2 lightflavors. The innermost 
characteristic line coincides with the prediction from lattice QCD [13] for the critical line at small 
fj,, which is represented by the dashed line with the hatched error band. 



The characteristic Une emanating from Tc is naturally related to the critical 
line Tc(/i) enclosing the hadronic phase. The comparison, in Figure 4, of our 
result for the curvature of the critical line at = 0, which can be calculated 
in lattice QCD [13], is a nontrivial and successful test of the extension of the 
quasiparticle approach to ^ > 0. The quark number susceptibility x{T) = 
dn/dfi\^=o is another quantity which has been computed on the lattice. As a 
second derivative of the pressure it is a very sensitive benchmark, and it agrees 
nicely with our result, see Figure 5. Finally, as a direct confirmation of our 
mapping procedure, the quasiparticle model can also successfully describe the 
available lattice data for p{fi,T) with 2 + 1 flavors [15]. It should be noted, 
however, that these direct calculations are so far restricted to small lattice sizes, 
resulting in still rather large uncertainties. With these supporting arguments we 
consider our results from the extrapolation of the lattice data at /i = (with 
controllable small uncertainties) as a realistic estimate for the equation of state 
at not too small temperatures. 

At low temperatures, matter will undergo a transition to a color-supercon- 
ducting state, with a different quasiparticle structure than presumed in our 
quasiparticle approach. Nonetheless, pairing affects the thermodynamic bulk 
properties only at the relative order of 0(A^//i^), where the estimated gap 
A < 100 MeV is comfortably smaller than the chemical potential. Therefore, 
our equation of state is a reasonable approximation even at small temperatures 
(maybe except for the pressure where it becomes very small). 
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Figure 5. The quark number susceptibility for Nf — 2, calculated from the quasiparticle 
model with the same parameters as in Figs. 3 and 4, for several chemical potentials compared to 
the lattice data [14] at = 0. 

Relevant for the following discussion is the equation of state in the form e{p), 
at T 0. Although both the pressure and the energy density deviate sizably 





Figure 6. The pressure and the energy density, scaled by the free results, for Nf — 2. 

from their ideal values, cf. Figure 6, we observe an almost linear relation 

e{p) =ap + 4:B , (12) 

as shown in Figure 7. The same linear form, with a = 3, follows in the 
familiar MIT-bag approach where p{x) = cx^ — Smit {"^ is /i or T, and c is 
some constant). Although this parameterization of the pressure is clearly in 
contradiction with the thermodynamic lattice results, the bag model relation for 
e{p) serves as a standard to compare our results with. We note in passing that 
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a linear relation between e and p, with a slope a « 3 is to be expected since 
both quantities scale with the forth power of either fj, or T. Causality requires 
furthermore a > 3. We find a 3.2, as a sign of the seen interaction effects, 
and B f» 3T^. 

3. Implications for quark stars 

Although there exist some lattice data for two light and one heavier quark, 
they still have larger uncertainties than for the case Nj = 2 considered so far. 
Nevertheless, based on the available data we can point out some interesting 
implications. From the universality of the scaled pressure p/pq{T/Tc) for 
various numbers of flavors, as observed on the lattice [16], we can expect that 
also at /i > the results do not change much besides a scaling from the different 
number of degrees of freedom in the physical case. On more general grounds we 
expect an approximately linear relation for e{p) for large ranges of the pressure. 

This expectation is confirmed by analyzing the available data, and approx- 
imating the lepton component in f3 equilibrated matter by an ideal gas. The 
results for e{p) ap + 4B, in particular, are found to be rather insensitive 
even under large arbitrary variations of the values A and Tg which parameterize 
g{p = 0, T). The slope parameter is found to be constrained by 3 < a < 3.5. 
The value of B, on the other hand, is in our approach directly linked to the tran- 
sition temperature, B 3T^^. Since Tc ~ 160 MeV is measured on the lattice, 
we have a definite prediction for the absolute scale in the equation of state. Our 
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estimate, B^^^ ?a 210 Me V, is substantially larger than the typical value of the 
bag constant obtained from fitting hadron spectra, -B^IT ~ 150 MeV. 

For a linear equation of state, the Tolman-Oppenheimer-Volkov equations 
imply a scaling property for the total mass and the radius of the star, 

M,R^B-^/^. (13) 



The effects of deviations of a from 3 being small, cf. Figure 8, we pointed 




Figure 8. The mass-radius relation, for tire linear equation of state (12), of a static quark star. 

out [11] the possibility of the existence of very dense and compact objects, 
M « O.DMsun and R Rr: 6 km, composed mainly of quark matter. Similar 
values for the maximal mass and radius were found in a perturbative approach 
with a physically motivated choice of the renormalization scale [17]. In a 
Schwinger-Dyson approach [18], M O.TMsun and R « 9km were obtained. 
It is interesting to compare these values with the mass and the radius of the object 
RX J1856: M ^ O.dMsnn and ii ?a 6km [19], which are not compatible with 
any hadronic equation of state. The precise values for M and R are, however, 
still under debate; for a recent discussion, we refer to [20]. 

Taking into account the effects of the hadronic crust of the star, its properties 
become sensitive to details of the hadronic equation of state, and of the transi- 
tion. For a strong first order transition, as suggested in [21], a new branch in 
the mass-radius diagram could exist, see Figure 9. The observation of so-called 
twins, as emphasized in [22], could be an exciting astrophysical indication of 
the existence of a non-hadronic phase in the center of neutron stars. 
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Figure 9. For a strong first order transition, a new stable branch can exist in the mass-radius 
relation (Figure from [17]). 
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